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We report on a study aimed to explore the dependence of the genuine multiparticle cor- 
^3 ■ relations on the number of sources when the influence of other possible factors during 

multihadron production is avoided. The analysis utilised the normalised cumulants cal- 
culated in three-dimensional phase space of the reaction e + e~— >Z° — > hadrons using a 
large Monte Carlo sample. The multi-sources reactions were simulated by overlaying a 
few independent single e + e~ annihilation events. It was found that as the number of 
sources S increases, the cumulants do not change significantly their structure, but those 
of an order q > 2 decrease fast in their magnitude. This reduction in the one-dimensional 
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rapidity cumulants can be understood in terms of combinatorial considerations of source 
mixing which dilutes the correlations by a factor of about 1/5" ?_1 . The diminishing of 
the genuine correlations is consistent with recent cumulant measurements in hadron and 
nucleus induced reactions and should also be relevant to other dynamical correlations like 
the Bose-Einstein one, in e + e~ — >W + W~ — > hadrons and in nucleus-nucleus reactions. 
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1 Introduction 



During the last decade an increase interest has been shown for the genuine multiparticle corre- 
lations in multihadron final states of hadronic, e + e~ and other reactions M. Recently OPAL, in 
its study of e + e~ annihilations on the Z° mass, has established the existence of strong genuine 
multihadron correlations up to the fifth order ||. In hadron-hadron, like proton-proton, colli- 
sions the correlations of more than three particles have also been observed . In contrast to 
this situation, in heavy ion collisions, at low energies and/or in reactions of light nuclei, genuine 
correlations are found to have non-zero values only up to the third order ||. Furthermore it 
has been found out that in general these correlations become weaker as the reaction average 
multiplicity increases. In nucleus-nucleus collisions at high energies, of tens and hundreds GeV 
per nucleon, the two-particle correlations are the only one that survive ||[7|-||]. 



This correlation dependence on the average multiplicity is very similar to the one observed 
in the investigations of multiparticle dynamical fluctuations, i.e. variation of many particle 
bunches in restricted phase space regions [Q. In these studies, known as intermittency analy- 
ses, the average multiplicity dependence has been proposed to be the consequence of a mixing 
of several independent emission sources [§, [U], [□]]. As a result, the dynamical fluctuations in 
nucleus- nucleus collisions are already well accounted for by two-particle correlations [f7|,^2|, 
whereas in hadron-hadron interactions P,IcSfl and in e + e~ annihilations || higher order genuine 
correlations do exist. 



An analogous situation is also observed in Bose-Einstein correlations (BEC) where identical 
bosons are correlated when they emerge from the interaction in nearby phase space. A genuine 
three-pion BEC has been detected in hadron-hadron reactions |14[ and found to be even more 
pronounced in e + e" 



annihilations 15 



In contrast to these, no genuine three-pion BEC were 
found in nucleus-nucleus collisions where the three-body correlations were well reproduced in 
terms of two-particle BEC Since the intermittency phenomenon and BEC seem to be 

closely related the dependence of many sources on the strength of the BEC cannot be 
excluded. The superposition of emitters may also be a reason for the suppression of BEC of 
hadrons produced from W-boson pairs in e + e~ annihilations at LEP2 energies where the over- 
lapping of hadrons affect the accuracy of the W mass measurements [17| . 



All this, as well as the obvious intrinsic interest in the genuine correlations which carry most 
of the dynamics of the hadron production process, points to the need of dedicated studies 
aimed to investigate the correlation dependence on the number of emission sources. Here we 
propose to study this dependence by grouping several e + e~ — > hadrons events to represent a 
multi-emission sources of particles. To obtain significant results, even when only few sources 
are considered, one needs a very high statistics, like that which can be supplied by a Monte 
Carlo (MC) generated sample, to minimise the calculation error and thus be sensitive to the 
correlations dependence on the number of sources. Another advantage of using MC generated 
events is in its possibility to generate a multihadron sample free from contamination of other 
processes like, for example, e + e _ — > t + t~ — > hadrons. 

In this Letter we present the results of a MC study on the effect of several emission sources 
on the genuine higher order multiparticle correlations. The study was based on a generated 
sample of about 5 x 10 6 events of the reaction e + e~ — > Z° — > hadrons which passed a full 
simulation of the OPAL detector at LEP and did reproduce rather well the measured genuine 



high order correlations present in the OPAL hadronic Z° decay data 0. Moreover, the e + e~ 
— > Z° — > hadrons annihilations should represent well the one emission source situation in 
contrast to events produced in hadron-hadron interactions. Our analysis on the dependence of 
many sources on the genuine correlations was thus carried out in a way that avoided effects of 
other reaction features, like the multiplicity which was discussed recently in connection with 
two-particle BEC analyses in [18| and [19]. In as much that final state interactions between 
hadrons coming from different sources can be neglected, our correlation study based on e + e _ — > 
Z° — > hadrons annihilations may be extended to other types of reactions since the hadronisation 
process is believed to rest on a common basis pi . 



2 The analysis method 



The analysis is based on a generated sample of hadronic Z° decays using the Jetset 7.4 MC 
program |21j including a full simulation of the OPAL detector at LEP [22|. The MC sample 



also included initial-state radiation and effects of finite lifetimes. The parameters of the pro- 
gram were tuned to yield a good description of the measured event shape and single particle 
distributions 



The selection criteria for multihadron events used here are identical to the ones previously 
utilised by OPAL in their recent data analysis of multiparticle correlations ||. In particular, 
selected events were required to have at least five charged tracks each having at least 20 mea- 
sured points in the jet chamber where the first point had to be closer than 40 cm from the 
beam axis. The cosine of the polar angle of the event sphericity axis with respect to the beam 
direction was required to be less than 0.7 to ensure the event to be within the volume of the 
detector. The sphericity axis was calculated by using all accepted tracks and electromagnetic 
and hadronic calorimeter clusters. 



To simulate several emission sources we did overlay several e + e~ — > Z° — > hadrons generated 
events and analysed the correlations between pions as if they were created in a single event. 
The kinematic variables are defined within each generated e + e~ event with respect to its own 
sphericity axis. For correlation analyses of variables like rapidity this procedure is equivalent to 
the one where the events are rotate to a common sphericity axis. This simulates multi-sources' 
events with hadrons, here all taken to be pions, emerging from a common emitter. To note is 
that in this procedure the average event multiplicity is directly proportional to the number of 
sources. In our analysis each generated event was used only once, and hence required a very 
large MC event statistics. 



To extract the genuine dynamical g-particle correlations, we used bin-averaged normalised 



factorial cumulant moments, or cumulants, first proposed in Ref. |24j as a tool for search for 
genuine multiparticle correlation, 
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where C q (yi, . . . , y q ) are the g-particle correlation functions given by the inclusive g-particle 
density distributions p q (yi, . . . ,y q ) in terms of cluster expansion, e.g., 
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£3(2/1,2/2,2/3) = 03(2/1,3/2,1/3) - X! Pi (2/1)^2 (2/2, 2/3) + 2 pi(yi)pi (2/2V1 (2/3) • (2) 

(3) 

Here M is the number of equal bins of a width by into which the event phase-space is divided 
and the subscript (3) denotes the number of permutations. For simplicity we show all formulae 
in one-dimensional (e.g., rapidity) phase space. 

The feature of the C g -functions is that they vanish whenever there are no genuine correlations, 
i.e., the correlations are due to those existing in lower orders. The correlations extracted are 
of a dynamical nature since the cumulants share with normalised factorial moments (the inter- 
mittency analysis tool) the property of statistical noise suppression. 

In this paper we computed the cumulants as they were used in experimental studies, in partic- 
ular we used the form applied in Ref. [@ namely, 

Kfi . Y M , fcM 

q T?L 1 N m (N m -i)-..(N m -q + l) ' K) 
Here, the factors are the unnormahsed factorial cumulant moments, or the Mueller mo- 



ments P5 |, calculated for the mth bin. These factors represent the correlation functions C q 
integrated over the bin and N m is the number of particles in the mth bin summed over all the 
M events. The definition (|3[) takes into account the non-uniform shape of the single-particle 
distribution and the bias when the cumulants are computed at small bins. 

The cumulant calculations were performed in the three-dimensional phase space of the kine- 
matic variables commonly utilised in this kind of studies [III], namely: 



• The rapidity, y = In J(E + p\\ ) /(E — p\\ ) , with E and p\\ being the energy and longitudinal 
momentum of the hadron in the interval —2.0 < y < 2.0; 

• The transverse momentum in the interval 0.09 < pr < 2.0 GeV/c; 

• The azimuthal angle, < $ < 2tt, calculated with respect to the eigenvector of the mo- 
mentum tensor having the smallest eigenvalue, in the plane perpendicular to the sphericity 
axis. 

These variables are defined with respect to the sphericity axis, in a way and within the intervals 
similar to those used in a recent OPAL analysis || and in other cumulant studies [1] . 



3 Genuine correlations and the number of sources 
3.1 Monte Carlo studies 

In Fig. 1 we compare the MC based cumulants of orders q = 2, 3 and 4 calculated from a single 
e + e~ — > Z° — >• hadrons source (solid symbols) with those obtained by overlaying seven such 
events representing seven hadronic sources (open symbols). The calculations were performed 
in one-dimensional rapidity, in two-dimensional rapidity vs. azimuthal angle subspaces and in 
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three-dimensional phase space of rapidity, azimuthal angle and transverse momentum. 
The following observations can clearly be made from Fig. 1. 

• The existence of a dynamical component, i.e. rise of the cumulants with increasing number 
of bins M, is seen to be present both in the single source as well as in the case of many 
sources. Although the slopes of this scaling behaviour are smaller for several sources than 
for a single source, they are still strongly present. It is also evident that the character 
of the scaling is kept as the number of sources increases e.g., for one source as well as 
for seven sources the one- dimensional (rapidity) cumulants level off at the same M value. 
No such saturation exists for the one and seven sources cumulants in the two and three 
dimensions. 

• The genuine dynamical correlations, measured by the cumulants, significantly decrease 
with the increase of the number of sources. This decrease is stronger for higher order 
correlations namely, whereas the two-particle cumulants suffer a reduction of an order 
of magnitude as the number of sources increases from one to seven, the four-particle 
cumulants diminish by three or four orders of magnitude. 

• The hierarchy of the K q cumulants is reversed as the number of sources increases. The 
cumulants derived from the single-source events increase with increasing g-order so that 
K% < K% < K4 , whereas the hierarchy in the cumulants calculated for seven sources 

(7) (7) (7) 

is reversed namely, K\ > K% > K\ . In addition, the multi-sources cumulants of order 
q > 2 have almost the same reduced value namely, w K± < 0(0.1). This last feature 
does not change as the dimension increases. 

• The overall dominant feature of the analysis results is the diminishing value of the higher 
order cumulants as the sources number increases leaving the K 2 to be the dominant 
genuine multiparticle correlations. 

The observed diminished correlations is further illustrated in Fig. 2 where the M-averaged 
one dimensional (rapidity) cumulants, are plotted with their errors against the number of 
sources. The cumulants were averaged over the M > 10 region where they are seen in Fig. 1 to 
approach an almost constant value. The values of these if* v are also listed in Table 1 together 
with the OPAL measured data cumulants of single e + e~ — > Z° — > hadrons events. These 
data cumulants for q < 3 are seen to agree with those derived from the MC sample. The one 
source MC based q = 4 cumulant lies lower than the measured data value but is still consistent 
within errors. 

It is obvious from Fig. 2 that the M-averaged rapidity cumulants of order q > 2 decrease fast 
with the increase of the number of sources. Already for two sources the hierarchy changes and 
the two-particle correlations visibly dominate over the higher order ones. At higher number of 
sources the dominant role of the two-particle correlations is even more pronounced. 

3.2 Correlation dilution due to source mixing 

Within the procedure adopted here for the simulation of multi-source events, it is clear that if 
a genuine correlation exists it can only be detected in groups of q pions emerging from the very 
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same source. In those g-group combinations which emerge from at least two sources, genuine 
correlations should not be present. This means that for K q cumulants that are calculated over 
all possible g-pion groups, the higher the number of sources the more diluted will be the signal 
for genuine correlations. 

For a given g-order the genuine correlation dilution factor is thus: 

Rq = (PG + P NG) > ( 4 ) 

where denotes the number of q— particle groups, e.g., pairs or triplets of pions, which emerge 
from the same source. The term P^ G stands for the number of all possible combinations of 
q— particle groups which emerge from at least two sources. 

Since all sources are produced in the same reaction and at the same energy, they do have an 
identical average charged multiplicity. For the estimation of R q we assume that all the S sources 
have the same fixed charged multiplicity n. In this case one has P 9 G = and the dilution 

factors at q = 2, 3 and 4 are given by 



R-2 



GO S n»l 1 



05 + 
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(4) S n»3 1 



(J) 3 + n* (J) + 3 @ (f ) + @ 2 (f ) + 2 n (j) (f ) S ' 



where the denominators include the number of all possible g-particle combinations in S sources 
of charged multiplicity n. This dilution factor dependence on the number of sources can also 



be derived in terms of cumulants ITC 



From these R q relations one can show that as long asn>g one obtains a general expression 
for the dilution factor, 

To compare the dilution factors R q with our correlation results shown in Fig. 2, they do have to 
be multiplied by K q v ^ which is a measure of the genuine g-order correlation present in a single 
sources. The solid lines shown in Fig. 2 thus represent the diluted cumulants = x R q . 

The striped areas in which the lines are embedded are the allowed regions when q is not ne- 
glected with respect to the multiplicity n. The agreement between the cumulant calculations 
and the dilution factors predictions is really remarkable for q = 2 and 3 and certainly is still 
well within the rather large errors of the q = 4 cumulants. 



For the order q = 2 one can relax the fixed charged multiplicity assumption and allow them 
to be different and still retain the R q ~ l/S^ 1 relation as long as the multiplicity distribution 
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is of a Poisson nature. This however is not the case for orders higher than 2. Nevertheless 
for order q = 3 the relation i? 3 = 1/S 2 , derived from the fixed multiplicity assumption, is 
still valid as it describes well the K^ v values up to at least thirteen sources (see Fig. 2). The 
large cumulants' errors associated with the q = 4 order prohibits to judge how accurate is the 
i?4 = 1/S 3 relation. 

An additional interesting and useful application of the relation R q ~ 1/S q ~ l is that it offers a 
method to estimate the average number of sources (S) via the cumulant averaged values over 
the large M region of two sequential g-orders through the ratio, 

( s ) - ^av(i) x jF^r ■ ( 7 ) 

3.3 Comparison with hadron and nucleus induced reactions 

As is already mentioned in the introduction, the genuine correlations measured in e + e~ an- 
nihilations are found to be weaker in hadronic interactions |3|-|5|, |13|| and even more so in 



nuclear collisions [@, In nucleus-nucleus collisions at ultra- relativistic energies only the 

second-order correlations were so far detected 

In Table 2 we list the results obtained by several experiments ||, §, [| on the M-averaged 
rapidity cumulant values for q = 2 and 3. These average values were taken over the M-region 
where the cumulants are seen in the published figures to reach a constant level. The hadronic 
reactions and their cumulants values are ordered according to their reported mean charged 
multiplicity, from the lowest value to the highest one. 

Table 2 shows that in hadron including nucleus induced reactions the two-particle correlations 
decrease rather fast as the mean multiplicity increases. However the three-particle correlations 
are found to be essentially non-existing even at moderately small mean multiplicity. Notwith- 
standing the possibility that production of hadrons in e + e~ annihilation may well be simpler 
than in hadron induced reactions, it may nevertheless be instructive to relate our findings 
to the measured correlation data listed in Table 2. In as much that the mean multiplicity 
increases with the number of sources, the decrease in the two-particle correlations and the 
absence of three-pion correlations in nucleus induced reactions is consistent with our findings 
which demonstrated the dilution of the correlations with increased number of sources. A quan- 
titative comparison between our findings and the correlations in nucleus-nucleus and hadron 
induced reactions is hard mainly because of the large errors associated with the average cumu- 
lant values. In particular the application of relation (|7]) is prohibited because most of the K^ v 
are consistent within errors with zero. 

Recently the two-pion BEC have been studied []BJ in pp reaction at centre of mass energy 
of 630 GeV as a function of multiplicity by using the normalised cumulants method similar 
to the one used here. In that analysis it has been found that the correlations of the cumu- 
lants of the like-sign pions as well as the opposite-sign pions decrease with the multiplicity n c h. 
From our analysis we expect the pair correlation to decrease as 1/S, where S is the number of 
sources. This indicates that indeed the multiplicity is at least partially proportional to S. The 
BEC dependence has also been investigated in the framework of the totally coherent emission 
picture |26[ and in the quantum optical approach where the conclusions were that these 
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correlations are weaker as the multiplicity increases. 



4 Summary and conclusions 

To investigate the effect of many emission sources on the genuine correlations in mutihadron 
final state we adopted a procedure which should minimise the confusion introduced by other 
variables like charged multiplicity. For the genuine correlations measurement we utilised the 
normalised cumulant method. To simulate the situation of many sources event we did overlay 
Monte Carlo generated hadronic Z° events treating them as one event. This Jetset7.4 MC 
sample of some five million events, tuned to the OPAL data taken at LEP1 on the Z° mass, has 
previously described rather well the measured correlations in the e + e~ — > Z° — > hadrons data. 

The results obtained here show that the cumulants, obtained from a single-source events and 
from events of many sources, almost do not change their structure with the decrease of the 
width of phase-space bins. This means that the scaling is preserved although larger slopes are 
seen to be in the case of one source compared with those for several sources. However, when 
the number of sources increases, the cumulants of order higher than two are suppressed and 
diminish to zero due to source mixing. The two-particle cumulants are also somewhat reduced 
in their value but they are way above the higher order ones and they are seen to completely 
dominate when the source number S exceeds the value ten. 



The one-dimensional (rapidity) correlations are very well reproduced by assuming that genuine 
correlation of the order q can only be present when all the q hadrons are emerging from the same 
source. Therefore the dilution of the genuine correlation signal is proportional to the ratio of 
the probability that the q hadrons will come from the same source. From simple combinatorial 
considerations this probability is approximately equal to l/S^ 1 . Thus a measurement of the 
one-dimensional correlation for two sequential orders renders the number of sources. 



The genuine correlations measured in hadron and nucleus induced reactions do follow qualita- 
tively the findings of our work. In particular in nucleus- nucleus reactions, where many sources 
are expected to contribute to the final hadronic state, the q > 2 orders are very small and 
indeed consistent with zero. The q = 2 order is still present but it is also getting smaller as 
the atomic number of the nuclei increases. The general decrease of the second order cumulants 
with the increase of multiplicity re-confirms the belief that the higher the multiplicity the larger 
the number of sources. Our results may also be useful for the understanding of other types of 
measured correlations like the Bose-Einstein interference of two and more identical bosons. It 
has been previously pointed out [57] that in the absence of final state interactions the BEC 
of the e + e~ — > W + W _ — > hadrons will be half of that of the Z° decay to hadrons. From our 
study it follows that the two-particle BEC, or any other correlations, in the two-source reaction 
e + e~ — > W + W _ — > hadrons should be reduced by a factor two as compared to that of the 
hadrons emerging from one W-boson. 
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Table 1: The Monte Carlo M averaged rapidity K^ v cumulants obtained in the e + e~ — > Z° 
— ► hadrons reaction compared with those obtained in a recent OPAL measurement of single 
data events. The averages were taken over the M region where the cumulants approached a 
constant value (see Fig. 1). 



No. of K™ (M > 10) 



sources <? = 2 q = 3 q = 4 



Sample 



0.45 ±0.01 0.67 ±0.04 1.36 ±0.21 Data g 

"1 0.486 ±0.002 0.632 ±0.017 0.950 ± 0.225 MC 

2 0.260 ±0.001 0.171 ±0.006 0.147 ±0.034 

4 0.124 ±0.001 0.041 ±0.003 0.016 ± 0.012 

7 0.072 ±0.001 0.013 ±0.003 0.003 ± 0.008 

13 0.037 ±0.001 0.004 ±0.003 0.001 ±0.009 



Table 2: The M averaged rapidity cumulants K^ v , of orders q = 2 and 3 measured in several 
hadronic reactions. The cumulant values were averaged over the M regions where they were 
seen to approach a constant value. These quoted values were estimated from the relevant 
published figures given in the references listed in the table. 



/ \ Kf d 4-- Beam energy D . 

(rich) = 2 5 Reaction toJ Ref. 

q = 2 g = 3 (GeV) 

"^8 0.32 ±0.02 0.26 ±0.12 vrp 250 §^ 

21.1 0.21 ±0.04 0.14 ±0.18 pEm 200 [|* 

> 50 0.34 ±0.02 0.12 ±0.03 AuEm 10.6 A 

73.3 0.21 ±0.03 0.05 ±0.07 SiEm U.5A [|,|] 

81.1 0.20 ±0.05 0.00 ±0.12 OEm 60A 

154.9 0.11 ±0.05 0.01 ±0.18 OEm 2QQA 

216.1 0.25 ±0.05 0.08 ±0.10 SEm 2QQA 

272.6 0.09 ±0.05 0.00 ±0.18 SEm 2QQA 

289.8 0.11 ±0.08 0.02 ±0.10 SAu 2QQA 
355.0 0.08 ±0.05 0.00 ± 0.08 SAu 200,4 

383.9 0.07 ±0.04 0.00 ± 0.05 SAu 200A 



* Calculations are based on the measured factorial moments. 
** Semicentral collisions. 
*** Central collisions. 
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Figure 1: The MC predicted cumulants of order q = 2, 3 and 4 as a function of M 1//D , where 
M is the number of bins of the D- dimensional subspaces of the phase space of rapidity (y), az- 
imuthal angle ($>), and transverse momentum (pr)- The solid symbols represent the cumulants 
for a single source, while the open symbols are the cumulants values of seven sources. 
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Figure 2: The dependence of the averaged rapidity cumulants K^ v of order q = 2, 3 and 4 
on the number of e + e~ sources. The cumulants were averaged over the M -range where they 
approached a constant value (see Fig. 1). The lines represent the expected dilution according 
to Eq. U^) where q is neglected in comparison to the multiplicity n. The striped areas are the 
allowed regions when q is not neglected with respect to n (see text). 
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